Europhysics Letters 



PREPRINT 



Osmotic Pressure of Solutions Containing Flexible Poly- 
mers Subject to an Annealed Molecular Weight Distribu- 
tion 

A. MlLCHEV 1 , J. P. WlTTMER 2 (*), P. VAN DER SCHOOT 3 and D. LANDAU 4 

1 Institute for Physical Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bul- 
garia 

2 Departement de Physique des Materiaux, Universite Claude Bernard and CNRS, 
69622 Villeurbanne Cedex, France 

3 Department of Applied Physics, Technische Universiteit Eindhoven, 
Postbus 513, 5600 MB Eindhoven, The Netherlands 

4 Department of Physics & Astronomy, University of Georgia, Athens, Ga 30602, USA 



PACS. 61.25.Hq - Macromolecular and polymer solutions. 
PACS. 82.35.+t - Polymer reaction and polymerization. 

PACS. 61.25Hq - Macromolecular and polymer solutions; polymer melts; swelling. 



Abstract. - The osmotic pressure P in equilibrium polymers (EP) in good solvent is investi- 
gated by means of a three dimensional off-lattice Monte Carlo simulation. Our results compare 
well with real space renormalisation group theory and the osmotic compressibility K oc (f>A(j)/AP 
from recent light scattering study of systems of long worm-like micelles. We confirm the scal- 
ing predictions for EP based on traditional physics of quenched monodisperse polymers in 
the dilute and semidilute limit. Specifically, we find P oc </f> 2 ' 3 and, hence, K oc <?r - 3 in the 
semidilute regime — in agreement with both theory and experiment. At higher concentrations 
where the semidilute blobs become too small and hard-core interactions and packing effects 
become dominant, a much stronger increase is evidenced and, consequently, the compressibility 
decreases much more rapidly with <f> than predicted from semidilute polymer theory, but again 
in agreement with experiment. 



Introduction. Superficially, it would seem that the colligative properties of quenched and 
annealed polymers must differ vastly, because the former have a fixed molecular weight distri- 
bution, whilst the latter are in equilibrium with each other, and continually exchange mate- 
rial Q. In this paper, where we study by means of computer simulations the osmotic pressure 
of "equilibrium polymers" (EP) in a good solvent, we demonstrate that systems of EP, as 
widely varying as giant polymer-like surfactant micelles, and supramolecular aggregates of 
dyes, dipolar colloids and proteins behave in essence like conventional polymers. The 
work presented here complements earlier simulation studies, focusing on the density distribu- 
tion of sizes in dilute, semidilute and concentrated solution, covering both cases with j^] and 
without ring closure [EL pi]. Our motivation for studying the osmotic pressure is that it is a 
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physical quantity that is more readily accessible experimentally than a distribution function or 
a mean size, in particular in the regimes where the chains strongly overlap. In addition, since 
conformation space renormalisation group (RG) predictions for the pressure are available j(| , 
as well as experimental data of the osmotic compressibility covering at least two concentration 
regimes , a fairly comprehensive comparison between theory, simulations and experiment is 
possible. The agreement we find is remarkably good, strongly suggesting that our simulations 
not only accurately describe the size distributions of EP, as in fact already shown in previous 
work , but also capture their colligative properties correctly. 

After very briefly presenting the simulation method, we first present our results for the 
pressure, and discuss these in the framework of the scaling theory of polymers. Both systems 
with and without ring closure have been investigated. (We recall that in the magnetic analogy, 
used to model EP, rings are suppressed B.) Our focus will mainly but not exclusively be on the 
findings for the EP with suppressed ring closure. The reason is that ring closure is thought 
not to play a significant role in giant micellar systems J2|, [?) , for which the most extensive 
experimental data are available. Also, for reasons to be discussed, although differences in the 
behavior of EP with and without rings do show up in the dilute regime, the scaling functions 
remain for all practical purposes the same. Also presented in this Letter is the scaling behavior 
of the mean length of the EP, on the one hand because it is directly related to the chemical 
potential of polymerizing materials and on the other because it is subject of controversy || . 
The combined pressure and chemical potential data fully describe the thermodynamics of 
the system in hand. We stress the success of the ideas borrowed from the scaling theory of 
polymers in the dilute and semidilute regimes, but also show that in the high-density limit 
the properties of EP are dominated by packing effects, and not by chain connectivity or 
self-assembly. 

Algorithm and some technical details. As in our previous studies [p[0], we use an off-lattice 
Monte Carlo approach, generalizing an efficient bead-spring model for polymers of fixed length. 
For the direct, non-bonding interactions between the beads we use the repulsive part of a Morse 
potential. For all practical purposes one may view these beads as hard spheres of diameter 
d ~ 0.8, and effective volume v = ird 3 /6 « 0.25. Here, and below, units of length are given in 
units of the maximum bond length of the bonding potential between two neighboring beads 
on a chain. This bonding interaction is described by the so-called FENE (finitely extensible 
nonlinear elastic) potential. In addition to the FENE potential we assign a scission energy 
E > to every bond, to be paid whenever a bond is broken [Q. This scission energy is 
presumed to be independent of the chain length and the density. In our presentation we 
choose energy units such that the Boltzmann factor ks — li and work at fixed temperature 
ksT = 1. The two operational parameters which determine the system properties are the 
number density of beads (j>, and the scission energy E. Both have been varied over a wide 
range, E from 3 to 13, and <j> from 2 -8 to 2 to produce a sufficient chain length and density 
variation allowing for a critical test of the theoretical predictions 0,0, |f|. We estimate that 
in our simulations the effective volume fractions V(f> probed range from 10 -3 to 0.45. The 
latter value has to be compared with the (only slightly larger) hard-sphere freezing volume 
fraction of about one half. Therefore, our simulations extend all the way up to the melt regime 
of a dense liquid. The combined FENE and Morse potentials fix the mean bond length at 
I w d, essentially independent of chain length and density |j| . In our model the chains may 
be regarded as fully flexible 0]. 

The osmotic pressure and its scaling with concentration. As is usual in off-lattice Monte 
Carlo simulations, the pressure P is obtained by evaluating the virial, and adding to that the 
contribution ksTcf) from the kinetic degrees of freedom Jn|. In Fig. |] we present the results 
for the bare pressure of the self-assembled polymers as a function of the number density <fi, 
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Fig. 1 - Pressure P versus the number density cj> for systems of monodisperse Morse beads (crosses), 
equilibrium polymers with ring closure suppressed (open symbols) and equilibrium polymers with ring 
closure allowed (full symbols). Results for various scission energies E are indicated by the symbol 
shapes listed in the graph. The dashed and full lines indicate the predicted scaling in the dilute and 
semidilute regimes respectively, and the dotted line the ideal gas behavior of the Morse beads at low 
densities. 

Fig. 2 - Master curve of the reduced pressure y = P (N) / (j> versus the reduced mean aggregation 
number x = (N) jg ((/>), based on the notions of the scaling theory of polymers. Same symbols as 
in Fig. [j]. The bottom curve gives all the measured data points, the top curve (shifted for clarity) 
includes only those data for which (j> < 1. The dashed line gives the ideal gas law, and the bold line 
the power law, expected for a semidilute solution. The data collapse is almost perfect when the data 
taken in the high density regime (beyond the semidilute range) are excluded. In the top curve we 
compare our simulations with the estimate y = 1 + 4.3a; (thin line), and the RG theory prediction 
from H (dash-dotted line). Both curves virtually coincide. 



for various scission energies E. For reasons of comparison results are also presented obtained 
for a gas of monomeric beads, interacting solely via the Morse potential ("Morse beads"). 

For the Morse beads (indicated by crosses), we find ideal gas behavior at densities <j> <C 1. 
At about V(f> w 0.1 hard-core interactions become important, and P rises steeply. This is 
to be compared with our results for EP without ring closure (open symbols), and with ring 
closure (full symbols). Not surprisingly, switching on the polymerization in the EP systems 
acts to lower the pressure relative to that of the gas of Morse beads. Focusing on the EP 
without rings first, we observe that in the limit of low densities and low scission energies, i.e., 
in the weak aggregation limit, the pressure is comparable to that of the Morse beads, as it 
should (see the data for E = 3). As the scission energy and, hence, the mean chain length 
increases, distinct power laws are observed in the dilute and semidilute regimes. The bold 
dashed line compares the pressure in the dilute limit with what one expects on the basis of 
ideal gas behavior P = <f>/ (N), using the scaling of the mean chain length (N) oc <fi 0A6 in the 
dilute regime derived below. At higher densities the pressure becomes independent of E. In 
the limit of sufficiently small densities, but large mean chain length, the measured pressures 
converge towards P oc </> 2 3 . This is in line with the expected scaling behavior of polymers 
(quenched or annealed) in semidilute solution. In semidilute solution the pressure is governed 
by the number of blobs or entanglements per unit volume, independent of molecular weight. 
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It is therefore also independent of any self-assembly process. 

Of course, the blob picture is only valid for sufficiently low densities and large chain lengths. 
Our simulation data confirm this: the larger the scission energy E, the better the semidilute 
scaling regime is discernible, in essence because the chains then overlap at lower densities. At 
densities above <p w 1, roughly when the packing effects make themselves felt in the gas of 
Morse beads, the pressure of EP rises more strongly than predicted by scaling theory, but still 
remains independent of E. This is the melt regime discussed in the introduction. 

In Fig. ^| we attempt to construct a master curve covering all our data for the EP, essentially 
inspired by standard scaling theory for conventional polymers [^| , and elucidated in more detail 
in the next section. The vertical axis gives the pressure divided by the ideal gas pressure 
<p/(N). Given on the horizontal axis is the ratio of the mean aggregation number (N) and 
the number of monomers per blob g ~ ^r 1 ^ 3v ~ 1 \ where v ~ 0.588 is the Flory exponent of 
a self-avoiding walk in three dimensions. The bottom curve includes results from all densities 
probed, while the top curve excludes data for densities <j>>\, i.e., those data outside of the 
dilute and semidilute regimes. Not surprisingly, the data collapse is only successful if the data 
of the melt regime are excluded. 

Included in both figures are also results obtained for EP in which ring closure was allowed, 
and polymer rings compete with linear chains for the available monomers. (The data are 
indicated by the filled symbols for a single value of E = 7.) In a previous study j| we 
found that in the dilute limit most of the monomers are assembled in trimer rings (the lower 
cutoff), and that the polymerization transition and the crossover to the semidilute regime 
roughly coincide. The pressure data are consistent with this, since we find at small densities 
P = <fi/3. Beyond the crossover to the semidilute regime almost all of the additionally added 
monomers are included in linear chains. The fact that a finite amount of monomers remains 
bound in closed loops in the semidilute regime appears to be irrelevant, for we observe within 
statistical accuracy the same pressure with and without rings. The data collapse perfectly in 
Fig. |^, although a huge amount of trimer rings are present at low and intermediate densities. 
Apparently, and perhaps surprisingly [ [TT| , the presence of loops does not significantly affect 
the excluded volume screening. 

Theoretical considerations. Our findings for EP without rings can be understood from the 
following Ansatz for the grand potential density 

oo 

n[c(N)] = < N ) (M W) — I — fiN + E + f N (N, 0)) - R{4>), (1) 

N=l 

The sum over monomer mass N describes the self-assembly in terms of an ideal mixing entropy, 
a Lagrange multiplier or chemical potential \i ^ that ensures conservation of mass, and a 
free energy penalty E + /n{N, (jy) associated with the chain ends ||. Of the latter, Jn(N, (f>) 
renormalizes the bare scission energy E due to differences in the micro-environment of the ends 
and central parts of the chains. This term is different in the dilute, semidilute and concentrated 
regimes ||. Applying the well-known statistical properties of self-avoiding walks, we infer 
that in the dilute limit /at decreases logarithmically with the degree of polymerization, i.e., 
Jn(N) = — (7 — 1) log(A) with 7 « 1.158 the susceptibility exponent of the n — > vector model 
in three spatial dimensions [|||. Entering the semidilute regime, one expects excluded volume 
effects to be screened out, and Jm should level off to /j\r(0) = (7— l)(3v — 1) _1 log(0) p. It is 
not known theoretically how Jn (0) behaves in the melt regime, so we focus first our analysis 
to the dilute and semidilute regimes. 

The remaining contribution R(<fi) accounts for excluded volume interactions which depend 
only on density <f>, and are not conjugate to c(N). Minimizing eq.(Q) gives in the limit of long 
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Fig. 3 - Scaling of the reduced mean chain length (JV)/ exp(J5/2.93), versus the reduced density 
4>/ exp(— E/3.8). Only data taken in the dilute and semidilute regimes are shown (<J> < 1). The 
dash-dotted line indicates the prediction from RG theory [[|. Inset: concentration dependence of the 
free energy penalty /jv in the "melt" regime, as obtained directly from the simulations (triangles), 
and from the estimate /jv = P/<j> (line). 

Fig. 4 - Compressibility K = <j>dP/d4> versus the density <j>, from simulation (dotted lines), light 
scattering measurements (asterisks) and from RG theory of EP (dash-dotted line) ||. The bold 
dashed line gives the power law K oc <J> 1,46 in the dilute regime, the bold line on the right the power 
law K tx <J>~ 0,31 expected for a semidilute solution. At higher densities, in the melt regime, the 
compressibility drops much more rapidly due to packing restrictions. 



chains an exponential equilibrium distribution function c eq (JV) oc exp(— I?— /jv(JV, 0)+//JV) S, 
and an osmotic pressure P = — f2[c eg ] = <fi/{N) + R{4>)- Following De Genes S, we assume 
a constant free energy density of ksT w 1 per blob, giving P(N)/(j> = F(x) = 1 + x with 
x (x (N)/g(4>). This simple estimate can fit remarkably well the observed master curve of 
Fig. g, and is barely distinguishable from the more elaborate RG treatment of ref. J(|, also 
shown in the figure. Note that the functional eq. (|l|) may be readily generalized to include 
the possibility of ring closure ||, or even chain branching. However, as long as the excluded 
volume term R(<f>) remains dominated by entanglements, the equation of state should be 
independent of chain architecture, in line with our observations. 

Scaling behavior of mean chain length and chemical potential. To help analyzing the 
concentration and scission energy dependence of the pressure we evaluated the mean chain 
length (JV) as a function of these parameters, the results of which are given in Fig. ||. From 
the equilibrium length density distribution given in the previous section, we easily find that 
(JV) cx <j> a exp(SE) with a = 5 = 1/(1 + 7) w 0.46 in the dilute, and a = (1 + (7 - l)/(3f - 
l))/2 w 0.6, S = 1/2 in the semidilute limit 0,0]. This also fixes the chemical potential of 
the polymerizing material: — 1/fj, = ~f(N) in dilute, and —1/fi = (JV) in semidilute solution, 
at least for sufficiently large mean aggregation number. (The chemical potential has also 
been obtained directly from the exponential cutoff of the distribution functions, confirming 
these relations and yielding a similar scaling plot.) In Fig. ||, both the mean aggregation 
number and the density are scaled by their values at the crossover from the dilute to the 
semidilute regimes, obtainable by a matching procedure described elsewhere j|. The slopes 
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giving the a exponents in both regimes are clearly visible. Also included in the figure is the 
RG result of ||. As expected, the agreement is reasonable but not perfect, because the theory 
is based on an approximation valid in the vicinity of four spatial dimensions. Extrapolated 
to three dimensions the growth exponents become a ~ 0.44 in the dilute, and a « 0.56 
in the semidilute regime, slightly underestimating the results from our simulations in both 
regimes. The astounding accuracy of the scaling function for the pressure calculated from the 
RG theory, plotted in Fig. [|, is due to the circumstance that the inherent inaccuracies of the 
theory are hidden in the scaled pressure and the scaled mean chain length. 

As in Fig. ^ we have in Fig. |^ only included data for densities <j) < 1 to avoid entering 
the melt regime, where the scaling and RG theories are known to be no longer accurate J5|. 
This is illustrated in the inset of Fig. |^, where we have plotted for three values of E the 
concentration dependence of the free energy penalty /at, which in the simulations may be 
measured either from the relation (N) 2 cx <fiexp(E+ /jv) or from fitting the length distribution 
c( N) , which turns out to remain purely exponential in the melt regime at variance with recent 
lattice simulations []1 2f| attempting to support a questionable modeling assumption proposed 
by Bouchaud et al. [|1 3f| to rationalize the observed Levy flight dynamical behavior in giant 
micelles. The figure shows that in this regime the free energy penalty increases more or 
less linearly with concentration, /at(</>) ^ 4>- We speculate that packing effects give rise to 
the free energy penalty, disfavoring chain ends because they disturb the local structure of the 
chain fluid. A crude estimate for fjq can be gotten in terms of the work required to create two 
new chain ends against the ambient pressure. In the inset of Fig. |^, we compare /jy with this 
estimate PSV, presuming that the additional free volume SV « l/<fi associated with free ends 
is of the order of the mean monomer volume. Although our argument is certainly simplistic 
and further study is warranted, the qualitative agreement does suggest that packing effects 
rule at such high densities - this may of course also be inferred from Fig. |l|, which hints at 
the convergence of the monomer and EP pressures in the melt regime. 

Comparison with experiment: Osmotic compressibility. In Fig. |i| we compare our simu- 
lation results with the isothermal osmotic compressibility, K , of a particular giant micellar 
system, measured by means of light scattering JjJ . We deduce K from the computed pressures 
using the thermodynamic identity K cx (f>d(f>/dP. In order to minimize scatter and obtain a 
numerically meaningful differentiation, we first fit fourth order polynomials to each data set 
of log P versus log </>, and then differentiate these polynomials. Because the dilute regime 
is very small for large E, we have in addition extrapolated these data to lower densities us- 
ing the well tested relationship for the mean chain length in the dilute limit (see left side of 
Fig. H). In other words, we insist on K cx <p 1+a cx <fi 1A6 in dilute solution for large E. The 
experimental data have been shifted in such a way to match our semidilute regime where 
K cx 0~( 2 ~ 3,y )/( 3,y-1 ). This is allowed, because (i) the experimental data for K are not ab- 
solute and given in arbitrary units only 0, and (ii) our volume fractions are only effective 
ones. Reasonable agreement was found over all regimes for the scission energy E — 13. From 
this we conclude that we are able to simulate coarse-grained EP comparable in size to real 
giant micelles. Also included is the result from the RG calculation ||, which fits well to both 
experiment and simulation in the dilute and semidilute regime. Apparently, the semidilute 
regime is rather small in both simulations and in the experiments of Buhler and co-workers Q . 

Discussion. We have investigated the (osmotic) pressure and compressibility of solutions 
of equilibrium polymers. Results from off-lattice Monte Carlo simulation, renormalisation 
group theory, and recent light scattering data of giant micelles have been compared. At low 
density and sufficiently high scission energy, i.e., for long enough chains, the results confirm 
unambiguously the theoretical predictions based on standard polymer theory Q . 

The universal scaling of P and K with respect to density <f> and mean chain-length (N) is 
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shown to be identical to that of quenched polymers, in agreement with theory || . Specifically, 
we find that in the semidilute regime the power law scaling of the pressure P cx <fi 2 ' 3 and 
compressibility K oc </> -0 3 , in agreement with experimental findings in many solutions of 
giant micelles but only for sufficiently flexible ones Jl4|, [l5|. Whether or not rings are 
present does not seem to alter the scaling of the pressure in semidilute EP. Despite that there 
is consensus about the growth law in the dilute regime, now and again workers find growth 
exponents that deviate from the expected value close to one-half, varying from a « O.f to 
f .5 [p|, p~5|— p~7|| . From our simulations we put forward that for these systems additional physics 
must be involved, such as due to ring closure, branching, the presence of charges, or the 
presence of more than one aggregating component. Clearly, the issue calls for a more detailed 
experimental investigation. 

In our simulations the blob picture breaks down at volume fractions exceeding, say, 0.1. 
The pressure then increases more rapidly than with the 2.3 power law, although it does 
stay below that of the monomeric Morse beads at equal concentration. As a consequence, 
the deduced compressibility decreases more strongly in this limit than with the —0.3 power 
expected from scaling theory. Our interpretation is that the EP then enter the concentrated 
or melt regime. The light scattering data of Buhlcr et al. also point at a crossover to the 
concentrated regime 0. In real systems the crossover to the concentrated regime is non- 
universal, and likely to depend on both the persistence length of the aggregates as well as 
their diameter [ fL8| . We intend to investigate by means of computer simulation the influence 
of chains flexibility on the colligative properties of EP in the near future. 

* * * 

We thank M.E. Cates and J.-L. Barrat for stimulating discussions. 
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